The role of spin noise in the detection of nanoscale ensembles of nuclear spins 
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When probing nuclear spins in materials on the nanometer scale, random fluctuations of the 
spin polarization will exceed the mean Boltzmann polarization for sample volumes below about 
(100 nm) 3 . In this work, we use magnetic resonance force microscopy to observe nuclear spin fluc- 
tuations in real time. We show how reproducible measurements of the polarization variance can be 
obtained by controlling the spin correlation time and rapidly sampling a large number of indepen- 
dent spin configurations. A protocol to periodically randomize the spin ensemble is demonstrated, 
allowing significant improvement in the signal-to-noise ratio for nanometer-scale magnetic resonance 
imaging. 
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Magnetic resonance signals detected in conventional 
magnetic resonance imaging (MRI) experiments origi- 
nate from the slight alignment of nuclear spins induced 
by an external magnetic field. This thermal equilib- 
rium ( "Boltzmann" ) polarization gives rise to a mean 
fractional spin polarization that is typically quite small, 
AN/N < 10~ 4 . For the large ensembles of nuclear spins 
detected in MRI experiments, usually N > 10 15 spins, 
the Boltzmann polarization is the dominant source of 
spin alignment. However, as new techniques, such as 
magnetic resonance force microscopy (MRFM) [1, 2], 
push detectable volumes below (lOOnm) 3 , another type 
of polarization becomes increasingly important: the "sta- 
tistical" polarization. 

Statistical polarization arises from the incomplete can- 
cellation of randomly oriented spins. The instantaneous 
polarization, i.e. the difference AN between spin-up and 
spin-down populations, can be either positive or nega- 
tive and will fluctuate on a time scale that depends on 
the random flip rate of the spins (for example, due to 
spin-lattice relaxation). For a random ensemble of spin- 
1/2 nuclei, it follows from the properties of the binomial 
distribution that the statistical fluctuations have vari- 
ance <j\ n = N[l — (AN/N) 2 }, where the overbar indi- 
cates mean value. In the limit of small mean polarization, 
which is representative of most experiments, the variance 
simplifies to ct\ n = N. The existence of statistical po- 
larization was pointed out by Bloch in his seminal paper 
on nuclear induction [3], and has been observed experi- 
mentally by a number of techniques, including supercon- 
ducting quantum interference devices [4], conventional 
magnetic resonance detection [5, 6, 7], optical techniques 
[8], and MRFM [2, 9]. 

As detection volumes enter the nanometer-scale 
regime, the standard deviation of the polarization fluc- 
tuations (taw can easily exceed the Boltzmann polariza- 
tion AN = NfiB/ksT, where [i is the magnetic moment 
of the spin, B is the polarizing magnetic field and T is 
the temperature [10]. The dominance of statistical po- 



larization, as defined by ctatv > AN, occurs for sam- 
ple volumes V < (p,B / 'k^T)~ 2 p^ 1 , where p^ = N/V is 
the spin number density. Assuming conditions represen- 
tative of high-field MRI microscopy of protons in water 
(B = 10 T, T= 295 K and p N = 7-10 28 m" 3 ), the volume 
corresponds to ~ (230 nm) 3 . For MRFM detection of 19 F 
nuclei in calcium fluoride, as considered in this paper 
(B = 2.9 T, T = 4.5 K and p N = 5 • 10 28 m" 3 ), the vol- 
ume for statistical polarization dominance is ~ (37nm) 3 . 

Given that statistical polarization is such a strong fea- 
ture of nanoscale nuclear spin detection, it is worthwhile 
to consider efficient methods to harness it for imaging 
applications. In order to generate a signal that is pro- 
portional to the spin density of the sample, it is natural 
to consider using the variance as the "signal" . If the 
polarization AiV is measured for n independent configu- 
rations of the spin ensemble, the sample variance s\ N is 
estimated as 
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(The symbol a 2 is used to denote the true or theoretical 
variance, while s 2 is the estimated or "sample" variance.) 

The variance estimate s 2 ^ is subject to some uncer- 
tainty since only a limited number of independent spin 
configurations can be sampled. Textbooks on statistics, 
e.g. Ref. [11], show that the standard error of 
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If this "spin noise" is the only noise present, the overall 
signal-to-noise ratio (SNR) for the variance determina- 
tion depends only on n, 
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This equation reveals a basic strategy for statistical spin 
detection: one should rapidly sample as many inde- 
pendent spin configurations as possible. Independent 
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FIG. 1: (a) Experimental setup. A (4/im) 3 CaF2 par- 
ticle attached to the end of a single-crystal Si cantilever 
(fc = 2.6 kHz, k c = 86AtN/m, Q = 18,500) is positioned 
80 nm above a FeCo-magnetic tip and placed in a static ex- 
ternal field Bo = 2.85 T (_Bn||z). 19 F spins are excited by an 
rf magnetic field Bi = 4mT (£?i||y), generated by passing 
current I(t) through a microwire situated directly under the 
tip [16]. All experiments are carried out at 4.2 K and in high 
vacuum, (b) Protocol used for spin detection and randomiza- 
tion. Repetitive rf frequency sweeps with a center frequency of 
114.8 MHz invert the spins adiabatically twice per cantilever 
cycle [17, 18]. For spin randomization, the rf field is inter- 
rupted at fixed intervals r r by a burst of n/2 pulses. Spin 
polarization is completely uncorrelated before and after ran- 
domization pulses, (c) Signal processing protocol as described 
in the text. 



spin configurations can be obtained by periodically re- 
randomizing the ensemble. 

The above analysis of the SNR represents an idealized 
case. In real experiments, the spin polarization is mea- 
sured via an intermediate quantity, such as the magnetic 
force, and the measurements are corrupted by noise. In 
our experiments, the measurement noise, including the 
cantilever thermal noise, is significant, and the spin sig- 
nals have a finite correlation time. These factors must be 
taken into account when analyzing the signal statistics. 

As shown in Fig. 1, we observe the nuclear spin po- 
larization using an ultra-sensitive cantilever to detect the 
attonewton magnetic force between spins in the sample 
and a nearby magnetic tip [1, 13]. The cantilever me- 
chanical resonance is driven by cyclic adiabatic spin in- 
versions induced by rf frequency sweeps, and by thermal 
noise. A fiber-optic interferometer monitors the result- 
ing cantilever motion and the cantilever oscillation sig- 
nal is synchronously detected by a two-channel lock-in 
amplifier. The phase of the lock-in amplifier is set so 
that the spin signal plus thermal noise appear in the in- 
phase (X) channel, while only thermal noise is present 



in the quadrature (Y) channel. The lock-in signals are 
digitized and low-pass filtered in software so as to con- 
trol the overall measurement bandwidth. Since the spin 
signal and the measurement noise are statistically inde- 
pendent, the variance of the in-phase channel is given 
by (j\ = ofpjn + o^ oisc . The quadrature channel mea- 
sures only noise, so <7y = c 2 oisc . The spin portion of the 
variance is then given by of pin = rj| — a\ . 

In the experiment, the in-phase and quadrature chan- 
nels are recorded for a time r to tai and the spin signal vari- 
ance s 2 pin is estimated from the measured sample vari- 
ances according to s 2 pin = s 2 x — Sy. We now seek an 
expression for the SNR associated with s 2 pin that is sim- 
ilar to Eq. (3), but also includes measurement noise. To 
obtain statistically independent spin configurations, the 
spins are periodically randomized at intervals separated 
by time r r using bursts of short (7r/2) rf pulses. The 
spin polarization remains essentially constant between 
randomizations if r r <C r m , where r m is the correlation 
time of the naturally fluctuating spin polarization. In or- 
der to keep the analysis simple we assume that the X and 
Y channels are filtered using convolution-type niters that 
average the signals in the intervals between the periodic 
randomizations. The actual experiment used first order 
low pass filters instead of the assumed averaging filters, 
resulting in a slight reduction of SNR performance [14]. 

The output of the averaging niters arc assumed to 
be sampled synchronously with the randomizations, re- 
sulting in the measurement sequences Xi and Yi, where 
i = 1, ...,n and n ~ T to tai/''~r is the number of indepen- 
dent spin configurations. From these samples, the vari- 
ances s\ and Sy are calculated as in equation (1). The 
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Note that in the absence of measurement noise, (4) is 
consistent with the previous result in (3). 

To finish the analysis, we need to know how cr 2 oise de- 
pends on the randomization interval r r . For a convolu- 
tion filter that averages the signal over the randomiza- 
tion interval r r , the power transfer function is |G(/)| 2 = 
sin 2 (7r/r r ) /7r 2 / 2 r r 2 , which has an equivalent noise band- 
width of A/ = l/2T r (single-sided with units of Hz). The 
noise variance is given by cr 2 oisc = S noisc Af = S , noiso /2T r , 
where Snoisc is the single-sided power spectral density of 
the measurement noise in each channel. We assume this 
spectral density to be independent of frequency within 
the bandwidth of the filter. This is the case in our exper- 
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FIG. 2: Cantilever tip oscillation amplitude as a function 
of time, (a) In-phase signal X(t) in the presence of the 19 F 
spins fluctuating naturally with a spin signal correlation time 
of T m w 3.5 s. (b) Quadrature signal V(f) due to thermal 
fluctuations of the cantilever. The observed thermal noise 
correlation time (r « 200 ms) is set by cantilever and filter 
bandwidths. (c) Same as (a), but spins are periodically re- 
randomized every r r = 351 ms. (d) Normalized autocorrela- 
tion function of spin signals obtained by analysis of an 8 min 
long record and after subtracting the measurement thermal 
noise autocorrelation function. 



iment where the cantilever thermal noise has a relatively 
broad spectrum since the cantilever is strongly damped 
by feedback. The SNR can thus be written as 
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From the above equation, the SNR is found to be 
maximized when the randomization period is r r = 
S'noisc/v / 2<7g pin , which is equivalent to choosing a random- 
ization rate and associated filter bandwidth such that 
°noisc = ""spin/v^- With this choice of r r and assuming 
that Ttotai > r r , we find 
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We have experimentally demonstrated the advantage 



of periodic spin randomization while measuring statisti- 
cal polarization of 19 F in a CaF2 single-crystal sample. 
In Fig. 2(a) we show a typical 60 second duration record 
of X(t) from the lock-in amplifier as the spins are cycli- 
cally inverted by rf frequency sweeps. This output con- 
tains contributions from both the spins and the thermal 
noise, while the quadrature channel signal Y(t), shown 
in Fig. 2(b), is predominantly just thermal noise. The X 
channel variance is 0.32 A , which corresponds to a force 
variance of approximately 450 aN 2 . The Y channel vari- 
ance is much smaller, about 0.055 A , which corresponds 
to the thermal force variance of 78 aN 2 . Based on an es- 
timated lateral field gradient of 10 6 T/m, the observed 
spin fluctuations in the X channel correspond to an rms 
statistical polarization of about 1500 19 F spins. 

The spin fluctuations in Fig. 2(a) can be seen by in- 
spection to have long correlation times, on the order of 
seconds. By calculating the autocorrelation function as- 
sociated with this and other similar time records, the 
autocorrelation was found to be well fit by an exponen- 
tial decay with correlation time of r m ss 3.5 s [Fig. 2(d)]. 
T m is closely related to the intrinsic rotating-frame spin 
lifetime T\ p and also depends on a number of extrinsic pa- 
rameters, for example the amplitude and modulation of 
the rf field B\ [12, 16]. If only a single one minute wave- 
form, such as in Fig. 2(a), is used to estimate the variance 
of X(t), the error is large, approximately 35%, since the 
number of independent samples n = rtotalAin = 17 is 
relatively small. 

The signal X{t) obtained with periodic spin random- 
ization is shown in Fig. 2(c), where the randomization is 
achieved by periodically interrupting the cyclic adiabatic 
frequency sweeps by a burst of 7r/2 pulses, as shown in 
Fig. 1(c). For this example we use the randomization 
period r r = 351ms, which corresponds to a 7r/2-pulse 
burst every 500th adiabatic inversion. The burst con- 
sists of 20 pulses of 2 [is duration with center frequencies 
uniformly distributed over the same frequency range as 
the adiabatic sweeps. Judging from the autocorrelation 
function of the randomized fluctuations, the spins are 
very effectively scrambled by this protocol [Fig. 2(d)]. 
As expected, the autocorrelation function falls linearly 
and reaches zero at r r . 

The periodic randomization allows us to make an im- 
proved estimate of the spin signal variance. Based on 
a randomization repetition period r r = 351 ms and the 
same total measurement time of 60 s, the number of in- 
dependent measurements is now Ttotai/^r = 171, which 
results in an uncertainty in the variance of about 11%, 
an improvement of more than a factor of three. 

To more concretely illustrate the impact of periodic 
randomization on the SNR, we show in Fig. 3(a) a one- 
dimensional lateral imaging scan over the CaF2 object. 
Each point represents the estimated spin signal variance 
from a one minute data record. Two overall signal max- 
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FIG. 3: One-dimensional image scans in ^-direction of the 
calcium fluoride sample. The points are the measured vari- 
ances of the spin force and the solid lines are guides for the 
eye. Scatter of the data points in the top trace is due to the 
uncertainty of the variance estimates. Much improved SNR 
is seen in the lower trace where the spins are periodically 
randomized in order to rapidly probe a large number of inde- 
pendent spin configurations. Lock-in filter time constants for 
the top and bottom traces are t c = 700 ms and t c — 30 ms, 
respectively. 



ima are visible that result from the two lobes of the imag- 
ing point spread function [15]. For the scan with natural 
spin randomization, the large uncertainty in the variance 
results in large scatter in the scan data over regions where 
the spin signal is highest. The periodic spin randomiza- 
tion, performed with r r = 113 ms, results in dramatic 
reduction of the data scatter, especially for the points 
with the largest spin signal [Fig. 3(b)]. The slightly 
larger noise in the regions where there is no spin signal is 
due to increased measurement noise that results from the 
larger measurement bandwidth needed to accommodate 
the shortened spin correlation time. 

As indicated by Eq. (5) and shown in Fig. 4, we can 
vary r r in order to optimize the SNR. Here we measure 
the spin fluctuations at a position of large spin signal 
(the 670-nm position of the lateral scan shown in Fig. 
3). To determine each SNR data point, 50 one-minute 



spin 



data records were acquired. The sample variance s 
from each record w as cal culated and the S NR was de- 
termined by SNR = s 2 • /s,2 , where s 2 ■„ is the mean 

•J Spill/ spin Spill 

of the sample variances and s s 2 is the standard devia- 

spin 

tion of the variances. As shorter randomization intervals 
are used, the SNR increases until it reaches a maximum 
around r r = 50 ms. The SNR then drops for the short- 
est r r value tested, as expected, since the large associated 
measurement bandwidth is admitting more measurement 
noise than is optimal. Obtaining results for even shorter 
T r was hampered by the finite response time of the can- 
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FIG. 4: Data points showing SNR for one minute measure- 
ments as a function of the randomization interval r r . The 
solid line represents the expected SNR based on Eq. (5) us- 
ing known experimental parameters (no free fit parameters). 
At the optimum r r = 52 ms, spin noise power equals \/2-times 
the measurement noise power. 



tilever (12 ms for the cantilever damped by negative feed- 
back). At the optimum point, the achieved SNR is 12, 
or about 6 times better than the SNR achieved without 
any forced spin randomization. 

Figure 4 also plots the theoretical SNR as predicted by 
Eq. (5) using the following parameters determined from 
the experiment: of pin = 900 aN 2 and 5 no i se = 66aN 2 /Hz. 
The theory predicts an optimal SNR of roughly 13 for 
a one-minute measurement time when r r = 52 ms, in 
good agreement with the experimental results. Overall, 
the SNR values obtained experimentally are just slightly 
lower than predicted, probably because the filtering that 
was performed did not use the ideal convolution filter as- 
sumed in the analysis, but rather a first order low pass 
filter of equal noise bandwidth (time constant t c = r r /2). 

In conclusion, we have shown that noise inherent in 
the measurement of small, statistically polarized spin en- 
sembles can be mitigated by rapid re-randomization of 
the spins. Beyond the fundamental interest in the con- 
trol and measurement of spin ensembles in the nanome- 
ter regime, the practical impact on nanoscale magnetic 
resonance imaging is significant: the demonstrated 6x 
improvement in SNR allows a 36 x increase of imaging 
speed. 
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